Globular Cluster System Research Articles

Origin of the correlation between stellar kinematics and globular cluster system richness in ultradiffuse galaxies

ABSTRACT Observational surveys have found that the dynamical masses of ultradiffuse galaxies (UDGs) correlate with the richness of their globular cluster (GC) system. This could be explained if GC-rich galaxies formed in more massive dark matter haloes. We use simulations of galaxies and their GC systems from the E-MOSAICS project to test whether the simulations reproduce such a trend. We find that GC-rich simulated galaxies in galaxy groups have enclosed masses that are consistent with the dynamical masses of observed GC-rich UDGs. However, simulated GC-poor galaxies in galaxy groups have higher enclosed masses than those observed. We argue that GC-poor UDGs with low stellar velocity dispersions are discs observed nearly face on, such that their true mass is underestimated by observations. Using the simulations, we show that galactic star formation conditions resulting in dispersion-supported stellar systems also leads to efficient GC formation. Conversely, conditions leading to rotationally supported discs lead to inefficient GC formation. This result may explain why early-type galaxies typically have richer GC systems than late-type galaxies. This is also supported by comparisons of stellar axis ratios and GC-specific frequencies in observed dwarf galaxy samples, which show GC-rich systems are consistent with being spheroidal, while GC-poor systems are consistent with being discs. Therefore, particularly for GC-poor galaxies, rotation should be included in dynamical mass measurements from stellar dynamics.

Discovery of Globular Cluster candidates in the dwarf irregular galaxy IC 2574 using HST/ACS imaging

Abstract We report the discovery of 23 globular cluster (GC) candidates around the relatively isolated dwarf galaxy IC 2574 within the Messier 81 (M81) group, at a distance of 3.86 Mpc. We use observations from the HST Advanced Camera for Surveys (ACS) to analyse the imaging in the F814W and F555W broadband filters. Our GC candidates have luminosities ranging from −5.9 ≥ MV ≥ −10.4 and half-light radii of 1.4 ≤ rh ≤ 11.5 pc. We find the total number of GCs (NGC) = 27 ± 5 after applying completeness corrections, which implies a specific frequency of SN = 4.0 ± 0.8, consistent with expectations based on its luminosity. The GC system appears to have a bimodal colour distribution, with 30 per cent of the GC candidates having redder colours. We also find 5 objects with extremely blue colours that could be young star clusters linked to an intense star formation episode that occurred in IC 2574 ∼1 Gyr ago. We make an independent measurement of the halo mass of IC 2574 from its kinematic data, which is rare for low mass galaxies, and find log M200 = 10.93 ± 0.08. We place the galaxy on the well-known GC system mass-halo mass relation and find that it agrees well with the observed near-linear relation. IC 2574 has a rich GC population for a dwarf galaxy, which includes an unusually bright ω Cen-like GC, making it an exciting nearby laboratory for probing the peculiar efficiency of forming massive GCs in dwarf galaxies.

Imposters among us: globular cluster kinematics and the halo mass of ultra-diffuse galaxies in clusters

ABSTRACT The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest haloes as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyse the catalogues of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass haloes with M200 ≲ 1011.2 M⊙. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive haloes is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intracluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than 10 GC tracers are particularly prone to severe uncertainties. Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive haloes suggested by GC kinematics around some UDGs.

Constraining the stellar populations of ultra-diffuse galaxies in the MATLAS survey using spectral energy distribution fitting

ABSTRACT We use spectral energy distribution fitting to place constraints on the stellar populations of 59 ultra-diffuse galaxies (UDGs) in the low-to-moderate density fields of the MATLAS survey. We use the routine prospector, coupled with archival data in the optical from the Dark Energy Camera Legacy Survey, and near- and mid-infrared imaging from the Wide-field Infrared Survey Explorer, to recover the stellar masses, ages, metallicities, and star formation time-scales of the UDGs. We find that a subsample of the UDGs lies within the scatter of the mass–metallicity relation (MZR) for local classical dwarfs. However, another subsample is more metal-poor, being consistent with the evolving MZR at high redshift. We investigate UDG positioning trends in the mass–metallicity plane as a function of surface brightness, effective radius, axis ratio, local volume density, mass-weighted age, star formation time-scale, globular cluster (GC) counts, and GC specific frequency. We find that our sample of UDGs can be separated into two main classes: Class A: comprised of UDGs with lower stellar masses, prolonged star formation histories (SFHs), more elongated, inhabiting less dense environments, hosting fewer GCs, younger, consistent with the classical dwarf MZR, and fainter. Class B: UDGs with higher stellar masses, rapid SFHs, rounder, inhabiting the densest of our probed environments, hosting on average the most numerous GC systems, older, consistent with the high-redshift MZR (i.e. consistent with early-quenching), and brighter. The combination of these properties suggests that UDGs of Class A are consistent with a ‘puffed-up dwarf’ formation scenario, while UDGs of Class B seem to be better explained by ‘failed galaxy’ scenarios.

The galactic acceleration scale is imprinted on globular cluster systems of early-type galaxies of most masses and on red and blue globular cluster subpopulations

Context. Globular clusters (GCs) carry information about the formation histories and gravitational fields of their host galaxies. It was found before that the radial profiles of the volume number density of GCs in GC systems (GCSs) follow broken power laws, while the breaks occur approximately at the a0 radii. These are the radii at which the gravitational fields of the galaxies equal the galactic acceleration scale a0 = 1.2 × 10−10 m s−2 known from the radial acceleration relation or the MOND theory of modified dynamics. Aims. Our main goals here are to explore whether the above results hold true for galaxies of a wider mass range and for the red and blue GC subpopulations. Methods. We exploited catalogs of photometric GC candidates in the Fornax galaxy cluster based on ground and space observations and a new catalog of spectroscopic GCs of NGC 1399, the central galaxy of the cluster. For every galaxy, we obtained the parameters of the broken power-law density by fitting the on-sky distribution of the GC candidates, while allowing for a constant density of contaminants. The logarithmic stellar masses of our galaxy sample span 8.0 − 11.4 M⊙. Results. All investigated GCSs with a sufficient number of members show broken power-law density profiles. This holds true for the total GC population and the blue and red subpopulations. The inner and outer slopes and the break radii agree well for the different GC populations. The break radii agree with the a0 radii typically within a factor of two for all GC color subpopulations. The outer slopes correlate better with the a0 radii than with the galactic stellar masses. The break radii of NGC 1399 vary in azimuth, such that they are greater toward and against the direction to NGC 1404, which tidally interacts with NGC 1399.

Total mass slopes and enclosed mass constrained by globular cluster system dynamics

We study the total-mass density profiles of early-type galaxies (ETGs: ellipticals and lenticulars) with globular clusters (GCs) as kinematic tracers. The goal of this work is to probe the total mass distribution, parameterised with a double power-law profile, by constraining the parameters of the profile with a flexible modelling approach. To that end, we leverage the extended spatial distribution of GCs from the SLUGGS survey (⟨RGC,  max⟩∼8 Re) in combination with discrete dynamical modelling. We use discrete Jeans anisotropic modelling in cylindrical coordinates to determine the velocity moments at the location of the GCs in our sample. Assuming a Gaussian line-of-sight velocity distribution (LOSVD) and a combination of informative and uninformative priors we use a Bayesian framework to determine the best-fit parameters of the total mass density profile and orbital properties of the GCs. We find that the choice of informative priors does not impact the enclosed mass and inner slope measurements. Additionally, the orbital properties (anisotropy and rotation of the dispersion-dominated GC systems) minimally impact the measurements of the inner slope and enclosed mass. A strong presence of dynamically-distinct subpopulations or low numbers of kinematic tracers can bias the results. Owing to the large spatial extent of the tracers our method is sensitive to the intrinsic inner slope of the total mass profile and we find ᾱ = −1.88 ± 0.01 for 12 galaxies with robust measurements. To compare our results with literature values we fit a single power-law profile to the resulting total mass density. In the radial range 0.1–4 Re our measured slope has a value of ⟨γtot⟩= − 2.22 ± 0.14 and is in good agreement with the literature. Due to the increased flexibility in our modelling approach, our measurements exhibit larger uncertainties, thereby limiting our ability to constrain the intrinsic scatter σγ.

Do ultra diffuse galaxies with rich globular clusters systems have overly massive haloes?

ABSTRACT Some Ultra Diffuse Galaxies (UDGs) appear to host exceptionally rich globular cluster (GC) systems compared to normal galaxies of the same stellar mass. After re-examining these claims, we focus on a small sample of UDGs from the literature that have both rich GC systems (NGC >20) and a measured galaxy velocity dispersion. We find that UDGs with more GCs have higher dynamical masses and that GC-rich UDGs are dark matter dominated within their half-light radii. We extrapolate these dynamical masses to derive total halo masses assuming cuspy and cored mass profiles. We find reasonable agreement between halo masses derived from GC numbers (assuming the GC number–halo mass relation) and from cored halo profiles. This suggests that GC-rich UDGs do not follow the standard stellar mass–halo mass relation, occupying overly massive cored haloes for their stellar mass. A similar process to that invoked for some Local Group dwarfs, of early quenching, may result in GC-rich UDGs that have failed to form the expected mass of stars in a given halo (and thus giving the appearance of overly an massive halo). Simulations that correctly reproduce the known properties of GC systems associated with UDGs are needed.

Investigation of globular cluster of irregular galaxies

The objects of study are irregular galaxies and the corresponding systems of globular clusters (GCS). The purpose of this article is to compile a list of irregular galaxies containing GCS and find a statistical relationship between the physical parameters. Additionally, let’s aim to develop a method for constructing a non-stationary theory of their origin. It is noteworthy that only a limited number of results from observational data on the study of GCS in irregular galaxies have been published before. It has been identified 23 instances of GCS around irregular galaxies through observations. Based on observational data, lists of irregular galaxies containing the GCS were compiled. Im-type irregular galaxies are found to contain the largest number of GCS. Primarily, GCS are observed in Im-, Sm-, and dIr-type galaxies. The average number of clusters is dominated by Sm-type galaxies. In addition, dwarf irregular galaxies contain few clusters, and the absolute magnitude of their parent galaxies is smaller.
 Empirical relationships between the parameters of GCS and their parent galaxies are established. A linear dependence of the logarithm of the number of globular clusters in the system on the absolute magnitude of the host galaxy has been found. The study demonstrates that as the number of globular clusters increases, the mass of the host galaxy shows a linear growth. Additionally, a potential approach to the GCS origin theory has been explored in irregular galaxies by analyzing oscillation modes against the background of the collapse model

Scaling relations for globular cluster systems in early-type galaxies – III. The inner flattening of the radial distributions

ABSTRACT The current properties of globular cluster systems (GCSs) in the nearby Universe are a consequence of an evolutionary history ruled by accretion and disruption processes. The analysis of scaling relations for GCSs is crucial for understanding the relevance of different physical processes involved in their origin and subsequent survival. In this work, radial profiles are calculated for a sample of GCSs associated with nearby massive galaxies. This is supplemented with results from our previous papers, achieving a total sample of 53 GCSs with radial profiles homogeneously calculated, and spanning more than two orders of magnitude in stellar mass of the host galaxy. Scaling relations involving different parameters of the GCSs are obtained for the entire sample, and their significance in terms of different scenarios is considered. Finally, the results are analysed in the literature context.

Globular clusters contribute to the nuclear star clusters and galaxy centre γ-ray excess, moderated by galaxy assembly history

ABSTRACT Two unresolved questions at galaxy centres, namely the formation of the nuclear star cluster (NSC) and the origin of the γ-ray excess in the Milky Way (MW) and Andromeda (M31), are both related to the formation and evolution of globular clusters (GCs). They migrate towards the galaxy centre due to dynamical friction, and get tidally disrupted to release the stellar mass content including millisecond pulsars (MSPs), which contribute to the NSC and γ-ray excess. In this study, we propose a semi-analytical model of GC formation and evolution that utilizes the Illustris cosmological simulation to accurately capture the formation epochs of GCs and simulate their subsequent evolution. Our analysis confirms that our GC properties at z = 0 are consistent with observations, and our model naturally explains the formation of a massive NSC in a galaxy similar to the MW and M31. We also find a remarkable similarity in our model prediction with the γ-ray excess signal in the MW. However, our predictions fall short by approximately an order of magnitude in M31, indicating distinct origins for the two γ-ray excesses. Meanwhile, we utilize the catalogue of Illustris haloes to investigate the influence of galaxy assembly history. We find that the earlier a galaxy is assembled, the heavier and spatially more concentrated its GC system behaves at z = 0. This results in a larger NSC mass and brighter γ-ray emission from deposited MSPs.

Measuring distances to galaxies with globular cluster velocity dispersions

ABSTRACT Accurate distances are key to obtaining intrinsic properties of astronomical objects such as luminosity or size. Globular clusters (GCs) follow a well-defined relation between their absolute magnitudes and internal stellar velocity dispersions (σ), offering an independent way to measure distances to their host galaxies via high-resolution spectroscopy. This is reminiscent of the ‘Faber–Jackson’ for elliptical galaxies. However, unlike galaxies, GCs have a very narrow range of mass-to-light ratios and simple star formation histories. Here, we show that the GC MV−log10(σ) relation is linear, whose slope is identical for the Milky Way and M31 GC systems. Based on this, we use 94 Milky Way GCs which have distances from GAIA parallaxes, or proper-motion dispersion profiles to derive a ‘GC velocity dispersion’ distance (GCVD) to M31, obtaining (m − M)0 = 24.51 ± 0.08 (d = 798 ± 28 kpc), in excellent agreement with independent measurements. Combining data for these two galaxies to create a fiducial relation using 296 GCs with high-quality measurements, we obtain a zero-point uncertainty (±0.06 mag) corresponding to a distance uncertainty of $\sim 3~{{\ \rm per\ cent}}$. We then use GCVD to obtain a distance to the giant elliptical galaxy NGC 5128 (Centaurus A), finding (m − M)0 = 27.95 ± 0.09 (d = 3.89 ± 0.16 Mpc). This is in excellent agreement with, and in some cases more precise than, literature estimates from the tip of the red giant branch or surface brightness fluctuations. We apply GCVD to Local Group galaxies with appropriate data and find good agreement with literature values even in cases with only one GC velocity dispersion measurement.

Catalogue of model star clusters in the Milky Way and M31 galaxies

ABSTRACT Detailed understanding of the formation and evolution of globular clusters (GCs) has been recently advanced through a combination of numerical simulations and analytical models. We employ a state-of-the-art model to create a comprehensive catalogue of simulated clusters in three Milky Way (MW) and three Andromeda (M31) analogue galaxies. Our catalogue aims to connect the chemical and kinematic properties of GCs to the assembly histories of their host galaxies. We apply the model to a selected sample of simulated galaxies that closely match the virial mass, circular velocity profile, and defining assembly events of the MW and M31. The resulting catalogue has been calibrated to successfully reproduce key characteristics of the observed GC systems, including total cluster mass, mass function, metallicity distribution, radial profile, and velocity dispersion. We find that clusters in M31 span a wider range of age and metallicity, relative to the MW, possibly due to M31’s recent major merger. Such a merger also heated up the in-situ GC population to higher orbital energy and introduced a large number of ex-situ clusters at large radii. Understanding the impacts of galaxy mergers and accretion on the GC populations is crucial for uncovering the galaxy assembly histories.

The influence of globular cluster evolution on the specific frequency in dwarf galaxies

ABSTRACT Dwarf galaxies are known to exhibit an unusual richness in numbers of globular clusters (GCs), property quantified by the specific frequency (SN), which is high for dwarf and giant elliptical galaxies, but with a minimum for intermediate-mass galaxies. In this work we study the role that GC evolution has in setting this trend, for which we use N-body simulations to evolve GCs in dwarf galaxies and quantify their disruption efficiency. We selected five individual dwarf galaxies from a high-resolution cosmological simulation, which includes GC formation and follow-up of their paths inside the host galaxy. Then, the tidal history of each GC is coupled to nbody6++gpu to produce N-body models that account for both, the interaction of GCs with their galactic environment and their internal dynamics. This results in a GC mass-loss parametrization to estimate dissolution times and mass-loss rates after a Hubble time. GC evolution is sensitive to the particular orbital histories within each galaxy, but the overall result is that the amount of mass that GC systems lose scales with the mass (and density) of the host galaxy, i.e. the GC mass-loss efficiency is lowest in low-mass dwarfs. After a 12 Gyr evolution all simulated GC systems retain an important fraction of their initial mass (up to 25 per cent), in agreement with the high GC to field star ratios observed in some dwarfs, and supports the scenario in which GC disruption mechanisms play an important role in shaping the GC specific frequency in dwarf galaxies.

Comparing Globular Cluster System Properties with Host Galaxy Environment* * Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. These observations are associated with programs 14219 and 15265.

We present Hubble Space Telescope photometry in the optical (F475X) and near-infrared (F110W) bands of the globular cluster (GC) systems of the inner halos of a sample of 15 massive elliptical galaxies. The targets are selected from the volume-limited MASSIVE survey, and chosen to sample a range of environments from sparsely populated groups to the brightest cluster galaxies in dense clusters. We also present a quantitative model of the relation between (F475X − F110W) color and cluster metallicity [M/H], using simulated GCs. Because much of the GC population in such galaxies is built up through accretion, the metallicity distribution of the GC systems might be expected to vary with galaxy environment. The photometry is used to create a completeness-corrected metallicity distribution for each galaxy in the sample, and to fit a double Gaussian curve to each histogram in order to model the two standard red and blue subpopulations. Finally, the properties of the GC metallicity distribution are correlated against galaxy environment. We find that almost no GCS properties and host galaxy environmental properties are correlated, with the exception of a weak but consistent correlation between blue fraction and nth-nearest-neighbor surface density. The results suggest that the systemic properties of the GC systems, at least in the inner to midhalo regions, are influenced more strongly by the local environment at early times, rather than by the environmental properties we see today.

The X-Ray Binary-star Cluster Connection in Late-type Galaxies

We conduct one of the largest systematic investigations of bright X-ray binaries (XRBs) in both young star clusters and ancient globular clusters (GCs) using a sample of six nearby star-forming galaxies. Combining complete CXO X-ray source catalogs with optical Physics at High Angular Resolution in Nearby Galaxies-Hubble Space Telescope cluster catalogs, we identify a population of 33 XRBs within or near their parent clusters. We find that GCs that host XRBs in spiral galaxies appear to be brighter, more compact, denser, and more massive than the general GC population. However, these XRB hosts do not appear to be preferentially redder or more metal-rich, pointing to a possible absence of the metallicity-boosted formation of low-mass X-ray binaries (LMXBs) that is observed in the GCs of older galaxies. We also find that a smaller fraction of LMXBs is found in spiral GC systems when compared with those in early-type galaxies: between 8% and 50%, or an average of 20% across galaxies in our sample. Although there is a non-negligible probability of a chance superposition between an XRB and an unrelated young cluster, we find that among clusters younger than 10 Myr, which most likely host high-mass XRBs, the fraction of clusters associated with an XRB increases at higher cluster masses and densities. The X-ray luminosity of XRBs appears to increase with the mass of the cluster host for clusters younger than ∼400 Myr, while the inverse relation is found for XRBs in GCs.

Investigating the M GCS–M h Relation in the Most Massive Galaxies

The relation between the total mass contained in the globular clusters of a galaxy and the mass of its dark-matter halo has been found observationally to be nearly linear over five decades of mass. However, the high-mass end of this relation is not well determined from previous data and shows large scatter. We analyze the globular cluster systems (GCSs) of a homogeneous sample of 11 brightest cluster galaxies (BCGs) through DOLPHOT photometry of their deep Hubble Space Telescope images in the F814W filter. We standardize the definition of M GCS, the total GCS mass, by using the GC total population within a limiting radius of 0.1R virial, while the dark-matter halo mass M h is determined from the weak-lensing calibration of M h versus M bary. When these 11 BCGs are added to the previously studied homogeneous catalog of Virgo member galaxies, a total value for η = M GCS/M h is found to be (3.0 ± 1.8internal) × 10−5, slightly higher than previous estimates but with much reduced uncertainty. Perhaps more importantly, the results suggest that the relation continues to have a near-linear shape at the highest galaxy masses, strongly reinforcing the conclusion that accreted GCs make a major contribution to the GC populations at high galaxy mass.

Formation of globular clusters in dwarf galaxies of the Local Group

ABSTRACT The existence of globular clusters (GCs) in a few satellite galaxies, and their absence in majority of dwarf galaxies, present a challenge for models attempting to understand the origins of GCs. In addition to GC presence appearing stochastic and difficult to describe with average trends, in the smallest satellite galaxies GCs contribute a substantial fraction of total stellar mass. We investigate the stochasticity and number of GCs in dwarf galaxies using an updated version of our model that links the formation of GCs to the growth of the host galaxy mass. We find that more than 50 per cent of dwarf galaxies with stellar mass $M_{\star }\lesssim 2\times 10^7\, \mathrm{M}_\odot$ do not host GCs, whereas dwarfs with $M_{\star }\sim 10^8\, \mathrm{M}_\odot$ almost always contain some GCs, with a median number ∼10 at z  = 0. These predictions are in agreement with the observations of the Local Volume dwarfs. We also confirm the near-linear GC system mass–halo mass relation down to $M_{\mathrm{h}}\simeq 10^8\, \mathrm{M}_\odot$ under the assumption that GC formation and evolution in galaxies of all mass can be described by the same physical model. A detailed case study of two model dwarfs that resemble the Fornax dwarf spheroidal galaxy shows that observational samples can be notably biased by incompleteness below detection limit and at large radii.

The large-scale structure of globular clusters in the NGC 1052 group

ABSTRACT Prompted by the many controversial claims involving the NGC 1052 group, including that it hosts two dark matter-free galaxies with overluminous and monochromatic globular cluster (GC) systems, here, we map out the large-scale structure (LSS) of GCs over the entire group. To recover the LSS, we use archival optical CFHT imaging data. We recover two GC density maps, one based on universal photometric properties of GCs from simple stellar population models, and one based on the properties of spectroscopically confirmed GCs in DF2 and DF4 (the two dwarf galaxies with overluminous GC populations). Both selection methods reveal overdensities around the massive galaxies in the group, as well as around NGC 1052 itself, that are coincident with the positions of previously identified stellar streams and tidal features. No intragroup GCs are found connecting these structures to any of the dwarf galaxies. We find, however, two other dwarfs in the group hosting GC systems. These include RCP32 with two GCs with ages equivalent to the GCs around NGC 1052 and DF9 with three GCs with ages similar to the GCs around DF2 and DF4. We conclude that the GC distribution in the group does not strongly support any formation scenario in particular. It favours, none the less, scenarios relying on galaxy–galaxy interactions and on the coeval formation of GCs around the DM-free dwarf galaxies. These may include the recently proposed bullet-dwarf formation, as well as high-redshift tidal dwarf galaxy models.

Great balls of FIRE II: The evolution and destruction of star clusters across cosmic time in a Milky Way-mass galaxy

ABSTRACTThe current generation of galaxy simulations can resolve individual giant molecular clouds, the progenitors of dense star clusters. But the evolutionary fate of these young massive clusters, and whether they can become the old globular clusters (GCs) observed in many galaxies, is determined by a complex interplay of internal dynamical processes and external galactic effects. We present the first star-by-star N-body models of massive (N ∼ 105–107) star clusters formed in a FIRE-2 MHD simulation of a Milky Way-mass galaxy, with the relevant initial conditions and tidal forces extracted from the cosmological simulation. We select 895 (∼30 per cent) of the YMCs with >6 × 104 M⊙ from Grudić et al. 2022 and integrate them to z = 0 using the cluster Monte Carlo code, CMC. This procedure predicts a MW-like system with 148 GCs, predominantly formed during the early, bursty mode of star formation. Our GCs are younger, less massive, and more core-collapsed than clusters in the Milky Way or M31. This results from the assembly history and age-metallicity relationship of the host galaxy: Younger clusters are preferentially born in stronger tidal fields and initially retain fewer stellar-mass black holes, causing them to lose mass faster and reach core collapse sooner than older GCs. Our results suggest that the masses and core/half-light radii of GCs are shaped not only by internal dynamical processes, but also by the specific evolutionary history of their host galaxies. These results emphasize that N-body studies with realistic stellar physics are crucial to understanding the evolution and present-day properties of GC systems.

A Photometric Survey of Globular Cluster Systems in Brightest Cluster Galaxies

Hubble Space Telescope imaging for 26 giant early-type galaxies, all drawn from the MAST archive, is used to carry out photometry of their surrounding globular cluster (GC) systems. Most of these targets are Brightest Cluster Galaxies and their distances range from 24–210 Mpc. The catalogs of photometry, completed with DOLPHOT, are publicly available. The GC color indices are converted to [Fe/H] through a combination of 12 Gyr single stellar population models and direct spectroscopic calibration of the fiducial color index (F475W–F850LP). All of the resulting metallicity distribution functions (MDFs) can be accurately matched by bimodal-Gaussian functions. The GC systems in all of the galaxies also exhibit shallow metallicity gradients with projected galactocentric distance that average . Several parameters of the MDFs including the means, dispersions, and blue/red fractions are summarized. Perhaps the most interesting new result is the trend of blue/red GC fraction with galaxy mass, which connects with predictions from recent simulations of GC formation within hierarchical assembly of large galaxies. The observed trend reveals two major transition stages: for low-mass galaxies, the metal-rich (red) GC fraction f(red) increases steadily with galaxy mass, until halo mass M h ≃ 3 × 1012 M ⊙. Above this point, more than half the metal-poor (blue) GCs come from accreted satellites and f(red) starts declining. But above a still higher transition point near M h ≃ 1014 M ⊙, the data hint that f(red) may start to increase again because the metal-rich GCs also become dominated by accreted systems.